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Abstract 

In this paper, we study the problem of finding a hypersurface family 
from a given spatial geodesic curve in R 4 . We obtain the parametric 
representation for a hypersurface family whose members have the same 
curve as a given geodesic curve. Using the Frenet frame of the given 
geodesic curve, we present the hypersurface as a linear combination of this 
frame and analyse the necessary and sufficient condition for that curve to 
be geodesic. We illustrate this method by presenting some examples. 
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1 Introduction 

Geodesic is a well-known notion in differential geometry. A geodesic on a 
surface can be defined in many equivalent ways. Geometrically, a geodesic on 
a surface is the shortest path joining any two points of that surface. Geodesies 
are curves in surfaces that play a role anologous that of straight lines in the 
plane. A straight line doesn't bend to left or right as we travel along it [1]. 

In recent years, there has been various researches on geodesies. Kumar 
et al., [2] presented a study on geodesic curves computed directly on NURBS 
surfaces and discrete geodesies computed on the equivalent tessellated surfaces. 
Wang et al., [3] studied the problem of constructing a family of surfaces from a 
given spatial geodesic curve and derived a parametric representation for a surface 
pencil whose members share the same geodesic curve as an isoparametric curve. 
A practical method was presented by Sanchez and Dorado, [4] to construct 
polynomial surfaces from a polynomial geodesic or a family of geodesies, by 
prescribing tangent ribbons. Sprynski et al., [5] dealt with reconstruction of 
numerical or real surfaces based on the knowledge of some geodesic curves on the 
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surface. Paluszny, [6] considered patches that contain any given 3D polynomial 
curve as a pregeodesic (i.e. geodesic up to rcparametrization) . Given two pairs 
of regular space curves n (u), r 3 (u) and r 2 (v), r± (v) that define a curvilinear 
rectangle, Farouki et al., [7] handled the problem of constructing a C 2 surface 
patch R (u, v) for which these four boundary curves correspond to geodesies of 
the surface. Farouki et al., [8] considered the problem of constructing polynomial 
or rational tensor-product Bezier patches bounded by given four polynomial 
or rational Bezier curves defining a curvilinear rectangle, such that they are 
geodesies of the constructed surface. 

However, while differential geometry of a parametric surface in IR 3 can 
be found in textbooks such as in Struik [9], Willmore [10], Stoker [11], do Carmo 
[12] , differential geometry of a parametric surface in R™ can be found in textbook 
such as in the contemporary literature on Geometric Modeling [13, 14]. Also, 
there is little literature on differential geometry of parametric surface family in 
R n [3, 15, 16, 17], but not in M 4 . Besides, there is an ascending interest on 
fourth dimension [18, 19, 20]. 

On the other hand, various visualization techniques about objects in 
Euclidean n-space (n > 4) are presented [21, 22, 23]. The fundamental step to 
visualize a 4D object is projecting first in to the 3-space and then into the plane. 
In many real world applications the problem of visualizing three-dimensional 
data, commonly referred to as scalar fields arouses. The graph of a function 
/ (x, y, z) : U C M 3 — > R , U is open, is a special type of parametric hypersurface 
with the parametrization (x,y, z, f (x,y, z)) in 4-space. There is an existing 
method for rendering such a 3-surface based on known methods for visualizing 
functions of two variables [24] . 

In this paper, we consider the four dimensional anologue of the problem 
of constructing a parametric representation of a surface family from a given 
spatial geodesic in Wang et al. [3], who derived the necessary and sufficient 
conditions on the marching-scale functions for which the curve C is an iso- 
geodesic, i.e., both a geodesic and a parameter curve, on a given surface. We 
express the hypersurface pencil paramctrically with the help of the Frcnet frame 
{T, N,Bi,B 2 } of the given curve. We find the necessary and sufficient con- 
straints on the marching-scale functions, namely, coefficients of Frenet vectors, 
so that both the geodesic and parametric requirements met. Finally, as an ap- 
plication of our method one example for each type of marching-scale functions 
is given. 

2 Preliminaries 

Let us first introduce some notations and definitions. Bold letters such as a, 
R will be used for vectors and vector functions. We assume that they are smooth 
enough so that all the (partial) derivatives given in the paper are meaningful. 
Let a : I C M — > R 4 be an arc-length curve. If {T,N,B 1 ,B 2 } is the moving 
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Frenet frame along a, then the Frenet formulas are given by 



T' = fciN, 

N' = -fciT + feBi, (1) 
B; = -fc 2 N + fc 3 B 2 , 
B 2 = -fc 3 Bi, 

where T,N,B 1; B 2 denote the tangent, principal normal, first binormal and 
second binormal vector fields, respectively, ki (i = 1, 2, 3) the ith curvature func- 
tions of the curve a [25]. 

From elementary differential geometry we have 

a'{s) = T(s), 

a" (a) = fci(s)N(s), ( 2 ) 
M«) = KMII- 
By using Frenet formulas one can obtain the followings 

a"'(s) = -kfT(s) + k' 1 N(s) + k 1 k 2 B 1 (s), 
a (m) (s) = -3k 1 k' 1 T(s) + (-kf + k'{-k 1 kf)N(s) (3) 
+ (2k[k 2 + k x k' 2 ) Bi (s) + hk 2 k 3 B 2 (s) . 

The unit vectors B 2 and Bi are given by 

= a' (a) to" (a) W (a) 
21 J \\a> (s) xa" (s) xa'" (s)W' [ ' 

Bi (s) = B 2 (s)iT(s)iN(s). 

Since the vectors T, N,B 1 ,B 2 are a right-handed orthonormal, the sec- 
ond curvature k 2 and the third curvature fc 3 can be obtained from Eqn. (3) 

as 



Bi (s) ■ a'" (s) 

Ms) = h{s} ' (5) 

- B 2 (£) • (£) 
3(j M«)M«) ' 

where '•' denotes the standard inner product. 

Let {ei, e 2 , e 3 , e 4 } be the standard basis for four-dimensional Euclidean 

4 

space R 4 . The vector product of the vectors u = ^2mei, 

»=i 

4 4 

v =^«iei,w =^Wiei is defined by 

*=1 i=l 
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uxvxw 



ei e 2 e 3 e 4 

Ui U 2 u 3 114 
Vl V 2 V3 Vi 
Wl W 2 W3 W4 



[26, 27]. 

If u, v and w are linearly independent then uxvxw is orthogonal to 
each of these vectors. 



3 Hypersurface family with a common geodesic 

A curve r (s) on a hypersurface P = P (s, t, q) C M 4 is called an isoparamet- 
ric curve if it is a parameter curve, that is, there exists a pair of parameters to 
and qo such that r (s) = P (s, to, qo)- Given a parametric curve r (s), it is called 
an isogeodesic of a hypersurface P if it is both a geodesic and an isoparametric 
curve on P. 

Let C : r = r (s), L\ < s < L 2 , be a C 3 curve, where s is the arc- length. 
To have a well-defined principal normal, assume that r" ^ 0, L\ < s < 

Let T (s) , N (s) , B 1 (s) , B 2 (s) be the tangent, principal normal, first 
binormal, second binormal, respectively; and let k\ (s) , k 2 (s) and k 3 (s) be the 
first, second and the third curvature, respectively. Since {T (s) , N (s) , B 1 (s) , B 2 (s)} 
is an orthogonal coordinate frame on r (s) the parametric hypersurface 

P{s,t,q) : [L 1 ,L 2 }x[T 1 ,T 2 }x[Q 1 ,Q 2 ]^M. i 
passing through r (s) can be defined as follows: 



P (*,*,</) 



r (s) + (u (s, t,q),v (s, t,q),w (s, t,q),x (s, t, q)) 



( T W \ 

N(s) 
V B 2 (s) j 



Li < s<L 2 ,T 1 <t<T 2 ,Q 1 <q<Q 2 , 



where u (s, t 7 q) , v (s, t, q) , w (s, q) and x (s, q) are all C 1 functions. These 
functions are called the marching scale functions. 

We try to find out the necessary and sufficient conditions for a hyper- 
surface P = P (s, t, q) having the curve C as an isogeodesic. 

First to satisfy the isoparametricity condition there should exist to G 
Pi, T 2 ] and q G [Qi, Q 2 ] such that P (s, t , qo) = r (s), L\ < s < L 2 , that is, 



u(s,t ,q ) = v(s,t ,q ) = w(s,t ,q ) = x(s,t ,g ) = 0, 
to G [Ti,r 2 ], ? G[Qi,Q 2 ], ii<s<i 2 - 



(7) 
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Secondly, the curve C is a geodesic on the hypersurface P (s, t, q) if and 

only if the principal normal N (s) of the curve and the normal n (s, to, qo) of the 
hypersurface P(s,t, q) are linearly dependent, that is, parallel along the curve 

C [28]. The normal n(s,io,3o) of the hypersurface can be obtained by taking 
the vector product of the partial derivatives and using the Frenet formula as 
follows 



(s,t,q) ( du(a,t,q) \ 

-q s = ( 1 + g; v(s,t,q)k 1 {s)\T(s) 

+ (u (a, t, q) h (s) + dw(S d ^ q) - w {a, t, q) k 2 ( 8 )j N (a) 

+ (v (s, t, q) k 2 (a) + dw( -Q s t,q) - x (a, t, q) k 3 (s)j B x (s) 

+ (w(a,t,q)k 3 (a)+ dx{S ds t ' q) ^ B 2 (a), 



9P(s,t,q) 

dt 



dt 

dw(s,t,q) 
dt 



dt 



and 



d p (s,t,q) = du(s,t,q) T | dv{s,t,q) N 
dq dq ' dq 

+ 5w^) B (() + ^ (() 
dq dq 

By using Eqn. (7) we have 

dP (5, t , g ) ^ dP (a, t , q ) ^ dP (s,t ,q ) 

ds dt dq 

4>i 0, h, qo) T (s) - 2 (s, t , qo) N (s) 
+03 ( s ' *o> 9o ) Bi (s) - 4 (s, t , <Zo) B 2 (s) , 

where 



n(s,t ,<?o) = 



5 



01 (Mo,<?o) 



dv(s,tq,q ) d-w(s,t ,qa) dx(s,ta,qa) 

ds ds ds 

dv(s,t B ,q ) dw(s,t ,q ) dx(s,t ,q ) 

dt dt dt 

dv(s,to,qg) dvr(s, tp, go) dx(s,t ,q ) 

dq dq dq 



dv / <9w <9x dw <9x 



dw ( dv dx dv <9x 



ds \ dt dq dq dt 

dx ( dv dw dv dw \ 
+ 9l \~dt~dq~ ~ ~dq"dt) ' 



ds \ dt dq dq dt 



2 (s,i o ,3o) = 



1 + 



du(s,tg,q ) dw(s,tg,qg) dx(s,tg,qg) 



dw(s,t ,qo) 
dt 

9w(s,tp,qo) 



dx(s,t ,qg) 
dx(s,t ,qg) 



ds 

du(s,t ,q ) 
dt 

du(s,t ,q ) 
dq 

du \ ( dw dx dw dx \ dw 

ds J \ dt dq dq dt J ds 



dx ( du dw du dw \ 

ds \ dt dq dq dt J 



du dx du dx 

dt dq dq dt 



3 (Mo, go) 



1 + 



du(s,t a ,qo) dv(s,t a ,qa) 9x(s,t ,g ) 



ds ds ds 

du(s,t ,q ) dv(s,tg,qg) dx(s,t ,qg) 

dt dt dt 

du(s,tq,qa) dv(s,to,qg) dx(s,t ,qg) 

dq dq dq 

( du \ f dv dx dv dx \ dv / du dx 
ds J \dt dq dq dt J ds \ dt dq 
dx ( du dv du dv \ 
+ T)s~ \~~dt Yq ~ ~dq~~dt) ' 



<9u dx 

dq dt 



04(Mo,5o) = 



1 + 



du(s,tg,q ) dv(s,tg,q ) dw(s,t ,q ) 

ds ds ds 

du(s,t ,qo) dv(s,ta,qa) dw(s,t ,q a ) 

dt dt dt 

du(s,t ,q ) dv(s,t ,qa) d-w(s,t ,qg) 

dq dq dq 



= 1 + 



+ 



du \ ( dv dw dv dw 
ds J \dt dq dq dt 
dw ( dudv dudv\ 

ds \ dt dq dq dt J 



dv 



du dw 

dt dq 



du dw 

dq dt 



So, n (s, to, qo) 1 1 N (s) if and only if 



0i (Mo, 9o) = 03(Mo,5o) = 04 (Mo, 9o) = 0, 2 (s,io,5o) t^O, (8) 
to e [T U T 2 ], qo€[Qi,Q 2 ], ii<s<i 2 - 
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Thus, any hypersurface denned by Eqn. (6) has the curve C as an isogeodesic 
if and only if 



u(s,t ,qo) = v(s,t ,qo) = w(s,to,qo)=x(s,t ,qo)=0, 
0i(s,*o,9o) = 03 (Mo, 9o) = </>4 (Mo, go) = 0, cf> 2 (s,to,q ) ^ 0, (9) 
^ e [Ti,T 2 ], <foe[Qi,Q 2 ], £i<*<£2- 

is satisfied. We call the set of hypersurfaces defined by Eqn. (6) and satis- 
fying Eqn. (9) an isogeodesic hypersurface family. 

To develope the method further and for simplification purposes, we analyse 
some types of marching-scale functions. 

3.1 Marching-scale functions of type I 

Let marching-scale functions be 

u(s,t,q) = l(s)V(t,q), 

v(s,t,q) = m(s)V(t,q), 

vr(s,t,q) = n(s)W(t,q), 

x(s,t,q) = p(s)X(t,q), 

Li < s < L 2 ,T 1 < t < T 2 ,Qi < Q < Qz, 

whereZ(s), m(s), n(s), p(s), U(t,q), V(t,q), W(t,q), Xf^JeC 1 
and I (s) ^ ^ m (s) , n (s) ^ ^ p (s) , Vs <E [L 1 ,L 2 ]. By using Eqn. (9) the 
necessary and sufficient condition for which the curve C is an isogeodesic on the 
hypersurface P (s, t, q) can be given as 

U(t ,q ) = V(t ,qo) = W(t ,qo) =X(t ,qo) = 0, 



aV(to,go)5X(t ,go) 


dV {t 


,q ) <9X (t , 


go) 


dt dq 


dq 


dt 




dV(t Q ,q Q )dW (t ,q ) 


dV (to, 


go) dW(t Q , 


go) 


dt dq 


dq 


dt 




dW(t ,q )dX(t (h qo) 


dW (t Q 


,q ) dX. (t , 


go) 


dt dq 


dq 


dt 





- o, 

= 0, (10) 

* o, 



t G [Ti,T 2 ], g e [Qi,Q 2 ], 
By Eqn. (10) we have 9v ^' 9o) = and = 0. 
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So, Eqn. (10) can be simplified to 

U(t ,<Zo) = V(* ,<Zo) = W(< ,<Zo) = X(t ,« ) =0 
dV(t Q ,q ) dV(t ,q ) 



dt dq 
dW (t , g ) dX (t , qo) dW (t 0) g ) dX (t 0) go) 



9t 9g <9g <9f 

to G [T U T 2 ], q G [Qi,g 2 ]- 



0. (11) 

7^0, 



3.2 Marching-scale functions of type II 



Let marching-scale functions be 



u(s,t,q) = l(s,t)U(q), 

v(s,t,q) = m(s,t)V(q), 

w(s,t,q) = n(s,t)W(q), 

x(s,t,q) = p(s,t)X(q), 



Li<s<L 2 ,T 1 <t<T 2 ,Q 1 <q< Q 2 , 

where l(s,t), m(s,t), n(s,t), p(s,t), U (q) , V («) , W(g), X(g) G 
C^Also let us choose V (q ) = = U (q ) = = 0. By using Eqn. 

(9), the curve C is an isogeodesic on the hypersurfacc P (s,t,q) if and only if 
the followings are satisfied 

n(s,t )W{q ) = p{s,t )X(q ) = 0, 



{s,t ) dX(q ) dW{q ) dp (Mo) , n , 1<v . 
~dt M P (S ' <0) n (S ' ~dq — * °' (12) 

to e [Ti,T 2 ] , qo G [Qi,Q 2 ] , ii < 5 < L 2 . 
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3.3 Marching-scale functions of type III 



Let marching-scale functions be 



u(s,t,q) = l(s,q)U(t), 

v{s,t,q) = m(s,q)V(t), 

w(s,t,q) = n(s,q)W{t), 

x(s,t,q) = p(s,q)X(t), 



L 1 <s<L 2 ,T 1 <t<T 2 ,Q l <q< Q 2 , 

where l(s,q), m(s,q), n(s,q), p(s,q), 11(1), V (t) , W (t) , X (t) G 
C 1 . Also let V (t ) = = U(t ) = ^± = O.By usin Eqn. (9) wc derive 

the necessary and sufficient condition for which the curve C is an isogeodesic on 
the hypersurface P (s, t, q) as 

n (s, q ) W (t ) = p (s, q ) X (t ) = 0, 



dW (t ) dp(s,q ) dn(s,q ) . . . dX (to) 

n {a, t ) — —X (to) ~^~W (t ) p (s, q Q ) + 0, (13) 

to G [T U T 2 ] , q G [QuQi] , Li < * < L 2 . 

4 Examples 

Example 1. Let r (s) = cos (s) , | sin (s) , ^s, , < s < 2ir, be a curve 
parametrized by arc-length. For this curve, 



T(s) = r'(s) = ^-isin(s),icos(s),i , 

N(s) = (- cos (s) , - sin (s) , 0, 0) , 

r'( 3 )xr"0Qxr'"( 3 ) / V6 _ Vl\ 
21 j ~ ||r'( S )xr"( S )xr"'( S )|| ^''3' 3 )' 

Bi(«) = B 2 xTxN= |-^ s in(s),^cos(s),-^,-^ 



Let us choose the marching-scale functions of type I, where 
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and 



l(s) = m (s) = n (s) = p (s) = 1 



XJ(t,q) = (t-t )(q-q ), V (t, g) = 0, W(t,g)=t-t , X(f,g) = g-g , 

to € [0, 1] , g € [0, 1] , < s < 2tt. 
So we have 



u(s,t,q) = (i - t ) (g - g ) , 

v(s,t,q) = 0, 

w(s,t,q) = t-t , 

x(s,t,g) = g-g . 

The hypersurface 



P (s, t, g) = r (s) + u (s, t, g) T (s) + v (s, t, q) N (s) + w (s, t, q) B x (s) + x (s, t, q) B 2 (s) 
1 1 

= cos ( s ) _ 2 ^ ~~ ^ ^ ~ 90 ^ Sin ^ ~~ "2" ^ ~~ Sin ^ ' 
1 1 

- sin (s) + - (t - t ) (<? - 9o ) cos (s) + — (f - t ) cos (s) , 
1 1 

2 s + 2 (* ~ iQ ) ^ - «>) - -g- (g - go) , 

a/2 V2 , . , . \/6 , . V3 . .. 

— s + — (*- to) (q -<io)- — (t- to) - — (<?- go)) 

< s < 2tt, < t < 1, < g < 1, t e [0, 1] , g S [0, 1] , is a member of 
the isogeodesic hypersurface family, since it satisfies Eqn. (11). 

By changing the parameters to and go we can adjust the position of the 
curve r (s) on the hypersurface. Let us choose to = 5 and g = 0. Now the curve 
r (s) is again an isogeodesic on the hypersurface P (s, t, g) and the equation of 
the hypersurface is 



P(s,t,g) = 



(I C os (s) -1^-1^ (q + V3) sin (s) , 
\ sin (*) + \ (t - ^ (q + V3) cos (a) , 



1 \/2 



V3 
6 



1\ V6 
*"2 + 1T^ 



\/2 , V2 / 1 
^ S + - l*-2 



V6 /. 1 



>/3 



10 



The projection of a hypcrsurface into 3-space generally yields a three-dimensional 
volume. If we fix each of the three parameters, one at a time, we obtain three 
distinct families of 2-spaces in 4-space. The projections of these 2-surfaces into 
3-space are surfaces in 3-space. Thus, they can be displayed by 3D rendering 
methods. 

So, if we (parallel) project the hypersurface P (s, t, q) into the w = 
subspace and fix q = i we obtain the surface 



s,t, - 
' '8 



8\/3/ 1' 

*- g J sm ( s )- 



1 . , , l + 8^3 / :i , 

2 sm ( s ) + § I * - 2 ) cos W ' 



1 



1 

16 



V3 
6 



t - 



1 



V6 
24 j 



< s < 2tt, < t < 1 
in 3-space illustrated in Fig. 1. 



Fig. 1. Projection of a member of the hypersurface family with marching-scale 
functions of type I and its isogeodesic. 

Example 2. Given the curve parametrized by arc-length r (s) = ^| sin (s) , ^ cos (s) , 0, -^s^ , < 
s < 2ir, it is easy to show that 
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T (s) = r' (a) - ^ cos (a) , ~ sin (a) , 0, ^ j 

N(s) = (-sin(s) , -cos(s) ,0,0) , 
B 2 (.s) = (0,0,-1,0), 

Bi(s) = I— cos(s), — — sin (s) , 0, - - I . 

Let us choose the marching-scale functions of type II, where 
n (s, t) = s + 1 + 1, p (s, i) = (s + !)(*- to) , 



and 



So, we get 



t/ (g) = V (q) = 0,W(q)=q-q o ,X(q) = l. 



u(s,t,q) = 0, 

v(s,t,q) = 0, 

w(s,t,q) = (s + t + q)(q-q ). 

x(s,t,g) = (s + l)(t-to). 



From Eqn. (12) the hypersurface 



1 \/3 

P(M,<?) = (-sin(s) + — (s + t+l)(q-q )cos(s), 
1 \/3 

- cos (s) - — (s + t + 1) (g - g ) sin (s) , 
-(a + l)(t-to), 

5 (* + *+!) («-*>)) 

< s < 27T, < t < 1, < q < 1, is a member of the hypersurface family 
having the curve r (s) as an isogeodesic. 

Setting to = | and go = yields the hypersurface 
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P(s,t,q) = (-sin(s) + — (s + t + 1) qcos(s) , 
icos (s) - ^ (s + t + l)? sin (s) , 

By (parallel) projecting the hypersurface P (s, t, g) into the u> = subspace 
and fixing q=55o wc S°t the surface 

1 Vs 

(- sin (s) + (s + t + 1) cos (s) , 

1 

2 cos ^ ~ 1000 ( s + * + ^ sin ( s ) ' 
-(- + 1) 

where, < s < 2ir, < t < 1 in 3-space demonstrated in Fig. 2. 




Fig. 2. Projection of a member of the hypersurface family with marching-scale 
functions of type II and its isogeodesic. 



13 



Example 3. Let r (s) = sin (s) , i cos (s) , 0, , < s < 27r, be an 
arc-length curve. One can easily show that 



T (s) = r' (s) = ^ cos (a) , -± sin (s) , 0, ^ J 

N(s) = (-sin(s),-cos(s),0,0), 

B 2 (.s) = (0,0,-1,0), 

x, , n /V3 , N \/3 . , , „ l\ 

Bi(5) = I— cos(s), — — sin (s) , 0, - - I . 

for this curve. 

If we choose the marching-scale functions of type III, where 



and 



then 



n (s, q) = sin (s (q - q )) , p (s, q) = sq 2 
U{t)=V (t) = 0, W (t) = 1, X(t)=t-t 



u(s,t,q) = 0, 

v(s,t,q) = 0, 

w(s,t,g) = sia(s(q-q Q )), 

x(s,t,q) = sq 2 (t-t ). 

Thus, from Eqn. (13) if we take g ^ then the curve r (s) is an isogeodesic 
on the hypersurface 



P(s,t,q) = (-sin(s) + — cos(s)sin(s(g-g )), 
1 \/3 

- cos (s) - — sin (s) sin (s (q - q )) , 
-sg 2 (i - t ) , 
^fS-\sm(s(q-q ))), 

where tt < s < 3?r, < t < 1, < g < 1. 

By taking to = 1 and go = 1 we have the following hypersurface: 
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P(M,<?) = (-sin(s) + — cos(s)sin(s(gr-l)), 

1 \/3 

2 cos ( s ) _ — sin ( s ) sin ( s (l ~ !)) ' 
- S<Z 2 (t - 1) , 

I sin ( s ( g _ 1))). 

Hence, if we (parallel) project the hypersurface P (s, t, q) into the z = 
subspace we get the surface 



1 i/3 
P(s,q) = (-sin(s) + — cos (s) sin (s (q - 1)) , 

1 \/3 

- cos (s) - — sin (s) sin (s (g - 1)) , 
7r < s < 3-7T, < q < 1, in 3-space shown in Fig. 3. 



Fig. 3 Projection of a member of the hypersurface family with marching-scale 
functions of type III and its isogeodesic. 
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5 Conclusion 



We have introduced a method for finding a hypersurface family passing 
through the same given geodesic as an isoparametric curve. The members of the 
hypersurface family is obtained by choosing suitable marching-scale functions. 
For a better analysis of the method we investigate three types of marching-scale 
functions. Also, by giving an example for each type the method is verified. 
Furthermore, with the help of the projecting methods a member of the family 
is visualized in 3-space with its isogeodesic. 

However, there is more work waiting to study. For 3-space, one possible 
alternative is to consider the realm of implicit surfaces F (x, y, z, t) = and try to 
find out the constraints for a given curve r (s) is an isogeodesic on F (a;, y, z, t) = 
0. Also, the anologue of the problem dealt in this paper may be considered 
for 2-surfaces in 4-space or another types of marching-scale functions may be 
investigated. 
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